Mon. Not. R. Astron. Soc. 000, [THg] (2008) Printed 28 March 2008 (MN IaT E X style file v2.2) 



Ha long term monitoring of the Be star (3 Cep Aa 



oo 
O 

o 

(N 



00 

S 1 
6 



> 
in 

(N 

o 

oo 
o 



G. Catanzaro* 

INAF - Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy 



Accepted 2008 March 27. Received 2008 March 27; in original form 2008 February 12 



ABSTRACT 

Recent papers published in the last years contributed to resolve the enigma on the 
hypothetical Be nature of the hot pulsating star f3 Cep. This star shows variable 
emission in the Ha line, typical for Be stars, but its projected rotational velocity is 
very much lower than the critical limit, contrary to what is expected for a typical 
Be star. The emission has been attributed to the secondary component of the (3 Cep 
spectroscopic binary system. 

In this paper, using both ours and archived spectra, we attempted to recover the 
Ha profile of the secondary component and to analyze its behavior with time for a 
long period. To accomplish this task, we first derived the atmospheric parameters of 
the primary: T e a = 24000 ± 250 K and log g = 3.91 ± 0.10, then we used these values to 
compute its synthetic Ha profile and finally we reconstructed the secondary's profile 
disentangling the observed one. 

The secondary's Ha profile shows the typical two peaks emission of a Be star with 
a strong variability. We analyzed also the behavior versus time of some line width 
parameters: equivalent width, V/R, FWHM, peaks separation and radial velocity of 
the central depression. 

Projected rotational velocity (wsini) of the secondary and the dimension of the 
equatorial surrounding disk have been estimated, too. 



Key words: stars: emission-line, Be 
spectroscopic 



stars: individual: /3 Cephei - (stars:) binaries: 



X 



1 INTRODUCTION 

The star (3 Cephei (HD 205021, V = 3.2) is well known to be 
the prototype of a class of hot pulsating stars. Several studies 
reported in the recent literature have been outlined a star os- 
cillating with at least five frequencies: fi = 5.2497104 c d~ , 
f 2 = 5.385 c d" 1 and f 3 = 4.920 c d" 1 (lAerts et al ll 19941) an d 
U = 5.083 c d" 1 and f B = 5.417 c d" 1 (|Telting et al.ll 19971 ). 

Further, the star (/3 Cep A) is actually a complicated 
multiple system. In fact, beside to have a visual compan- 
ion (/3 Cep B, V = 7.9) at a distance of 13.4", it is also 
a member of a spectroscopic system whose second star 
(/3 Cep Aa) was dis covered using speckle interferometry by 
ICezari et alj (|l972l ) at a distance of «0.25". The parame- 
ter s of the close binary orbit have been determined later 
by IPigulskv fc Boratvrj (Il992l ) combining speckle interfer- 
ometry and the variations in the pulsation period, due to 
th e so-called light t ime e ffect. Recent speckle measurements 
bv lHartkopf et all (|200ll ) place the pos ition of the compan - 
ion to about 0.1" from the primary and lBalega et al.1 (2002) 
found a magnitude difference of 1.8 between the two com- 
ponents in a red filter centered on 810 nm and 60 nm wide. 
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More observations are still necessary to improve the orbital 
solution. 

Another characteristic that gives more interest to this 
object i s the v ariable emission observed in the Ha. Since 
iKarpovl (1 19331 ). the presence of this emission had b een 
reporting by seve r al au thors: IWilson fc Seddonl l|l956h or 
iKaper fc Mathiasi (|l995l ), who observed a decreasing emis- 
sion from 1990 to 1995. In the same years, f3 Cep has 
been observed with the coude spectrograph of the 2.6 me- 
ter telescope of th e Crim ean Astrophysical Observatory by 
IPan'ko fc Tarasovl (|l997t) . Over the entire period covered 
by their observations, the Ha was in emission with a pro- 
nounced two component structure and a small intensity 
above the continuum level. The same authors stated that in 
November 1987 the line showed an absorption profile with- 
out any sign of emission. 

The kind of variable emission discussed so far is typ- 
ically the characteristic of Be stars. Be stars are rapidly 
rotating B-type stars that lose mass in an equatorial, cir- 
cumst ellar disk (see the recent review bv lPorter fc Riviniusl 
(2003)). A general characteristic of Be stars is that they ro- 
tate very fast, typically at about 70%-80% of their critical 
limit (v cr it = \ j GM* /FL) corre sponding to a several hun- 
dred of km s" - 1 (|SlettebakHl982l ). On the contrary, the pro- 
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jected rotational velocity of fl Cep A (26 km s~ 1 . lMorel et al.l 
(2006)) is very much lower than this limit and this questions 
its r eal nature. 

lHadrava &: Harmaned (1 19961 ) separated emission and 
absorption components of the Ha profile from their 1996 
spectra and found that they move in antiphase. They spec- 
ulated on two possible scenarios: i) the observed emission 
arises from reprocessing of the light originating in the stel- 
lar photosphere in an outer envelope, it) the star is not a 
pulsating object but a rapidly rotating oblique rotator with 
a magnetic di pole geometry. 

Recently. ISchnerr et al.l (|2006l ) argued for the first time 
that the Ha emission observed from the /3 Cep system is not 
related to the slowly rotating primary star, but to the sec- 
ondary, being the latter a classical Be star. Unfortunately, 
they obtained only one spectrum at the Nordic Optical Tele- 
scope and, then, no conclusion about the variability of the 
emission could be drawn. 

In this paper we try to recover the recent emission his- 
tory of /3 Cep Aa. To achieve this goal, we firstly disentan- 
gled the spectra of the two components from the observed 
one. This has been done using a synthetic profile for the Ha 
of [5 Cep A, for which T e g and logy have been estimated as 
described in Sect. [3] Then, we analyzed the behavior with 
time of some Ha width parameters like its equivalent width, 
the full width at half maximum, the peaks separation and 
the radial velocity of the central depression. Finally, we at- 
tempted to estimate the v sin i of j3 Cep Aa and the dimen- 
sion of its surrounding disk. 



2 OBSERVATION AND DATA REDUCTION 

We observed f3 Cep for more than one year from July 2005 
to November 2006. To extend our coverage to the past years, 
we queried various astronomical archives on the internet. At 
the end, we collected a number of 72 spectra covering the 
period from September 1993 to November 2006, for a total of 
more than 13 years. In the following we detailed each single 
observatory and instrument: 

• 42 spectra have been taken from the archive of the Flit- 
ter Observatory, University of Toledo, OH USA, equipped 
with 1-m Ritchey-Chretien reflector. The echelle spectro- 
graph is connected to the Cassegrain focus of the telescope 
by fiber optic cable of 200 /im diameter. The detector is 
a CCD manufactured by EEV with 1200 x 800 pixels of 
22.5 /jm size. The spectra used in this paper have been ob- 
served with R ~ 26000; 

• 4 spectra have been downloaded from the Isaac Newton 
Telescope archive, three of them have been acquired with 
IDS spectrograph and one with MUSICOS. The lines of the 
Th-Ar lamp show that for all the spectra the resolution is 
Rw 33000; 

• 7 spectra have been downloaded from the archive of 
the Elodie@OHP 1.93 m telescope. Spectra have been re- 
duced by the stand ard pipeline procedure, described in 
Baranne et al.l (| 19961 ). Resolving power of these spectra is 
« 42000; 

• the 91 cm telescope of the INAF- Osservatorio As- 
trofisico di Catania, has been used by us to carried out 19 
spectra of our target. The telescope is fiber linked to a RE- 
OSC echelle spectrograph, which allows to obtain R = 20000 



spectra in the range 4300-6800 A. The resolving power has 
been checked using emission lines of the Th-Ar calibration 
lamp. Spectra were recorded on a thinned, back-illuminated 
(SITE) CCD with 1024 x 1024 pixels of 24 /im size, typi- 
cal readout noise of 6.5 e~ and gain of 2.5 e~/ADU. These 
spectra ha ve been extracted from a more complete set of ob- 
servations (|Catanzaro et all [200 8V in such a way to choose 
the one with the best SNR for each observing night. 

The stellar spectra, calibrated in wavelength and with 
the continuum normalized to a unity level, were obtained 
using standard data reduction procedures for spectroscopic 
observations within the NOAO/IRAF package, that is: bias 
frame subtraction, trimming, scattered light correction, flat- 
fielding, fitting traces and orders extraction and, finally, 
wavelength calibration. IRAF package rvcorrect has been 
used to include the velocity correction due to the Earth's 
motions, all the spectra were then reduced into the helio- 
centric rest of frame. 



3 ATMOSPHERIC PARAMETERS OF THE 
PRIMARY 

The first task we had to asses in our study was to compute 
the synthetic Ha profile of the primary component that we 
used in the reconstruction process of the secondary's profile. 
This calculation is possible if one knows the effective tem- 
perature and surface gravity of the primary component of 
the system. 

The approach we used in this study to determine T e g and 
loggr was to compare the observed and theoretical profiles 
of H7 and H/3 lines by minimizing the goodness-of-fit 
parameter: 

n z^y si obB j 

where TV is the total number of points, I b s and 7th arc, 
respectively, the intensities - normalized to the continuum 
- of the observed and computed profiles and <5/ b s is the 
photon noise. Errors in T c ff and logg are estimated as the 
variation in the parameters which increases the \ 2 of a unit. 

The observed profiles presented in Fig. [1] have been cal- 
culated by averaging a number of 2 5 spectra collected by u s 
during the Julian Day 2 554 041.45 (jCatanzaro et al.ll2008h . 
Doppler correction has been applied to remove any radial 
velocity shift due to pulsations. This average spectrum has 
been used in order to improve the signal-to-noise ratio of 
each single exposures, in particular we obtained SNR « 250 
in the continuum close to the H7 and SNR ~ 300 in the con- 
tinuum close to the H/3. 

Th eoretical profiles have b een computed with SYN- 
THEFI llKurucz fc Avrettl Il98ll ) on the basis of ATLAS9 
l|Kurucz I Il993l ) atmosphere models. To reduce the num- 
ber of free parameters, the rotational velocity of f3 Cep A 
has been determined fitting the line profiles of the Sim 
triplet at A 4552, 4567 and 4574 A. The value obtained of 
vsini = 30 ± 1 km s _1 is almost coincident with the line 

1 All the Kurucz codes (ATLAS9 an d SYNTHE) hav e been used 
in the Linux version implemented bv lSbordone et al.l d2004ri 
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Figure 1. Comparison between observed (black) and synthetic 
(red) profile of H7 and H/3 lines. Atmospheric parameters used for 
this calculation are: T off = 24000 ± 250 K and log g = 3.91 ± 0.10. 
For each Balmer line we reported also the residuals (O - C) ex- 
pressed in percentage. Note the strong oxygen lines present in 
the wings of the H7 line. 



Figure 2. In this figure we illustrated the method used to re- 
cover the Ha profile of the secondary. With bold solid line we 
represented the observed composite spectrum, with red dashed 
line we represented the synthetic primary component and with 
the normal line we indicated the reconstructed secondary profile. 



broadening of 29 km s _1 measured by I Morel et al.l (|2006l ) in 
their analysis of high resolution spectra (R = 50 000) . 

The best fit was achieved for the following adopted 
parameters: T cff = 24000 ± 250 K and log g = 3.91 ± 0.10. In 
Fig. [TJwe showed the comparison between the observed and 
computed lines profiles together with their residuals (O - C). 
It is interesting to note the two symmetrical bumps that ap- 
pear both on the blue and red side of H/3, approximately at 
± 150 km s _1 from the core of the line. They are located 
at the same velocity of the peaks observed in the Ha of 
/3 Cep Aa (see Fig. [2}, then they could be connected to the 
H/3 of the secondary component. On the other side, j3 Cep A 
is a non-radial pulsator, thus the peaks could be the result 
of the average process we undertook to recover the mean 
profiles used for T ott - and logg determination. This distor- 
tion, whatever its origin may be, could led to a maximum 
error on the parameters of the order of 2 %. 

In the particular case of our target, this procedure could 
be considered valid under the initial hypothesis that H7 
and H/3 are not influenced by the secondary component of 
the binary system. As a check on the consistency of our 
T c ff and log g, we estimated the abundance of silicon from 
lines of different ionization states. In particular we used 
the Sim triplet and the Siiv A4654 A. The abundance has 
been computed by means of spect ral synthesis analysis us - 
ing the microturbulence found by iGies fc Lambert! (|l992h . 
that is: 7.6 ± 1.5 km s~\ We derived -4.67 ± 0.08 from Sim 
and —4.65 ±0.10 from Siiv. The consistency between these 
two determinations convinced ourselves about the accuracy 
of our atmospheric parameters. As adopted Si abundance 
we computed the weighted average between the two previ- 



ous values, that is log N(Si) /N to t =-4.66 ± 0.06. This abun- 
dance, that could be expressed as log e{Si) = 7. 38± 0.0(0, 
is fully consistent with typical abundances for B-type field 
stars reported in various literature study: for instance, 
iDufton et all (|l99fj| ) gave loge(ffi) = 7.5 ± 02 for both the 
QB as sociations h and x P er an d Cep OB3. lGies fc Lamberd 
( 1992) found log e(Si) = 7.58 from their sample of B stars 
and lCunha fc Lambert (|l994l ) found 7.40 ± 0.15 in a sample 
of B stars in the Orion association. 

For the sake of comparison, we searched for other de- 
termina tions of T e tf and l ogff i n the literature. In par- 
ticular, iHevnderickx et al.1 (I1994T) found T eff = 24550 and 
log g = 3.772 and Niecmzura fc Daaszv'nska-Daszkiewwiczl 
(|2005h derived T e fi = 24150 ± 350 K and log g = 3.69. As the 
reader can see, despite effective temperatures are perfectly 
consistent with our value, our gravity is « 0.2 dex greater 
than previous determinations. Anyway, since other authors 
did not report their errors on loggr, no conclusion could be 
drawn about the consistency between our and literature val- 



4 EXTRACTION OF THE SECONDARY 
SPECTRUM 

The procedure we undertook to recover the Ha profiles of 
the secondary component could be summarized as follows: 



2 This conversion has been calculated considering normal helium 
abundance. 
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Figure 3. Ha profiles of /3 Cep Aa reconstructed with the procedure described in the text. In each box we reported the heliocentric 
julian day of the observation and a label indicating the observatory: Rit - Ritter Observatory, INT - Isaac Newton Telescope, Elo - Elodie 
and OAC - Catania Astrophysical Observatory. 



• first, we computed synthetic Ha profiles for the primary 
component using the atmospheric parameters derived in the 
previous section. Even for this line we used ATLAS9 and 
SYNTHE; 

• then, we recovered the secondary Ha line profiles from 
the following formula: 

-fprim T T F scc , 1 ^ 

i*tot = ; (1) 

1 + r 

where F to t is the observed profile, -F pr i m is the synthetic Ha 



profile of the primary, r is the luminosity ratio between the 
two compo nents computed fr om the magnitude difference 
A M — 1.8 (Ba lega et al.l l2002) opportunely scaled at the Ha 
wavelength. Then from Eq. |T| it follows: 

F scc = (1 + r) F tot - F plim r (2) 

Since the primary component is a pulsating s t ar wit h 
an amplitude of a b out 30 km s -1 jAerts et al.l (|l994h , 
iTelting et all l|l997T ). ICatanzaro et all (120081 ) ), it has been 
necessary to correct the synthetic spectrum for the radial 
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Figure 4. As in Fig. [3] 



velocity corresponding to the pulsating phase. For this aim 
we used the nearby Cn doublet at AA 6578, 6582 A, both 
belonging to f3 Cep A. For the sake of clarity, we show in 
Fig. [5] an example of our procedure. 

The profiles of f3 Cep Aa, recovered as described before 
and converted in the heliocentric velocity scale, are showed 
in Fig. [3] and in Fig. [4] Inspection of those two figures re- 
vealed the strong variability of /3 Cep Aa, even though no 
transition between shell into normal B spectrum has been 
observed. It always shows a double-peak symmetrical profile, 
with V/R ~ 1. Spectra are ordered by Julian Date. 



5 Ha VARIABILITY AND DISK PROPERTIES 

In this section we analyzed the behavior with time of some 
line width parameters of the Ha emission profile, i.e.: equiva- 
lent width (EW), the full width at half-maximum (FWHM), 
the peaks separation (AV), the radial velocity of the cen- 
tral depression (RVcd) and the ratio of blue to red peak 
intensities (V/R) computed after subtraction of the contin- 
uum levejfl In order to better measure positions and flux 



3 Since our spectra have been normalized to the unity, we com- 
puted this ratio as (V-1)/(R-1). 
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Figure 5. Variability of some quantities characteristics of the 
Ha emission profile plotted in function of the Heliocentric Julian 
Day. From bottom to top panel: equivalent width expressed in 
A, V/R, peaks separation AV and radial velocities of the cen- 
tral depression, all expressed in km s . Meaning of the symbols 
is: squares - data from Ritter Observatory, triangles - data from 
INT, asterisks - data from Elodie and circles - data from Catania 
Astrophysical Observatory. In the EW panel, the Ho emission in- 
crease going up; the dotted line in the V/R panel represents the 
V = R symmetry. The top axis is labeled in fraction of years. 



levels of blue and red peaks, we fitted each profile with a 
function defined as the sum of the two gaussians centered 
respectively on the blue and the red peak. A similar fit has 
been considered to compute the radial velocity of the central 
depression. All this data with their relative errors have been 
reported in Tab. IA1I 
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Figure 6. Correlation between V/R, expressed in logarithm, and 
radial velocity of the central depression. 



In Fig. [5] we plotted all these quantities versus time. As 
a general consideration, we observed two different regimes 
of variations: a smoothed one between Julian Date 2449250 
and 2453000 and a more rapid variation starting from 
JD > 2453000. 

In the bottom panel of Fig. [5]we reported the equivalent 
width of Ha versus time. By definition, it is straightforward 
that lines with emission profiles will have negative EW. In 
the plot we inverted the directions of y axis in such a way to 
have ascending curve in correspondence of greater emission. 
EW decreases from September 1993 to May 1999, then a 
new strong emission phenomenon occurs in July 2002 that 
remains constant till August 2003. Again EW decreases, it 
reaches a minimum during September 2005 and it increases 
again until November 2006, date of last spectrum of our 
series. 

As concerns the FWHM we found its value almost con- 
stant in time, about 400 ± 30 km s _1 , at least at the resolv- 
ing power of our spectra. 

The peaks separation reveled an increase from Septem- 
ber 1993 to May 1999, then it decreased reaching a minimum 
in June 2005. After that date, the separation begins to in- 
crease again (third panel from bottom of Fig. [5}. 

The fourth panel shows the behavior of the radial ve- 
locity of the central depression. It shows a progressive red 
shift till May 1999. Then, we observed a blue shift till June 
2005. After this date velocities present much more scatter, 
probably due to a worse quality of data. In general, from the 
variations of the radial velocity, one could draw some useful 
indication about the state of the disk: expanding, contract- 
ing or stationary. This is true provided that the radial ve- 
locity of the underlying star is well known. Unfortunately, 
since this is not possible with our data, we could not say 
anything about the equilibrium state of the envelope. 
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A similar trend has been observed for the V/R ratio, but 
with an exceptional increase (V/R > 1.5) for three julian 
dates: 51144.261 (Elodie), 51155.338 (INT) and 51325.682 
(INT). Cyclic variability o f V/R has already been observed 
in a number of Be stars (|Hanuschik et alIll995T ) and this 
evidence has been theoretically explained in terms of one- 
armed density waves proce ssing around th e central star, the 
so-called Okazaki's model (|Okazakilll99lh . When the dense 
part is on the approaching side of the disk, the V peak is 
higher than R. On the contrary, when dense part is in the 
receding part of the disk, R peak will be enhanced. 

It is we l l kno wn from literature (see for example 
iHubert et all (| 19871 )). that in addition to the V/R variation 
Be stars are always accompained by profile shifts: blueward 
when red component is stronger (V/R<1), redward when 
blue component is stronger (V/R> 1). This feature is again 
explained in term of the one-armed density model as pro- 
posed bv lOkazakil (|l996h . In Fig. we show the correlation 
between RV and V/R. It is clear that when V/R > 1 the pro- 
file is redshifted (RV > 0) and, on the contrary, the profile is 
blueshifted (RV< 0) when V/R< 1. 

For all these quantities, we attempted to search for pe- 
riodic variability u sing the Phase Dispersion Method (PDM) 
|StellingwerJll978l ) as coded in the NOAO/IRAF package, 
but no succesfully resu lts h ave been achieved . 

iDachs et al" | l|l986h and lHanusc hik (1986), studying the 
kinematics of the disks surrounding Be stars, found some 
useful correlations between line width parameters (in par- 
ticular, FWHM and peak separation) and stellar projected 
rotational velocity, v sin i. Later, t he existence of th o se cor - 
relations have been confirmed by IHanuschik et all (QUI). 
Using those correlations, we attem pted to derive the v sin i 
of the underlying star. In particular. IHanuschik et al.l (|l988T ) 
gave equations that linked each other FWHM, peaks separa- 
tion and vsini for objects whit EW of the Ha emitting disk 
region less than 15 A. From our measurements, we derived 
a velocity equal to ~ 230 km s _1 . 

The rotation law in Be star circumstellar disk is usu- 
ally written as u(r<j) oc r^ J , where the index j = 0.5 if 
for keplerian rotation and j= 1 for conservation of angular 
momentum. According to former author, the observed nor- 
malized peaks separation is a function of the outer radius of 
the disk, expressed in term of stellar radius, via the formula: 

~ T-* (3) 

In their paper, IHanuschik et al.l fl988) concluded that 
for their sample of B e stars the value of j is closer t o 1 than to 
0.5. On the contrary. iHummel fe Vranckenl ((2000) concluded 
that the formation of Ha shell profiles in Be stars generally 
requires that jssO.5, i.e. the keplerian rotation is a valid 
approximation. As for our object, according to the Okazaki's 
model, the correlation between RV and V/R showed in Fig. [6] 
is a strong argument for the validity of keplerian rotation. 
Thus, using Eq. [3] we estimated the radius of the outer disk 
surrounding j3 Cep Aa to be 3.38 R, (j = 0.5). 



6 CONCLUSION 

In this paper we analyzed the last 13 years of the Ha emis- 
sion history of /3 Cep Aa. Spectroscopic data, for a total of 



72 spectra, have been obtained both from INAF - Osserva- 
torio Astrofisico di Catania and from archives of other three 
observatories as described in Sec. [5] 

We used H/3 and H7 profiles to derive effective tempera- 
ture and surface gravity of the primary. By spectral synthesis 
we obtained: T cff = 24000 ± 250 K and log g = 3.91 ± 0.10. 

Regarding the secondary, we conclude that (3 Cep Aa 
is a Be star which shows a high degree of variability in Ha 
emission intensity as well as in the other line width param- 
eters here considered. From 1993 to 1999 we observed a de- 
crease in the emission and a subsequent increase till July 
2002, that remains constant till August 2003. After this long 
cycle extended for 10 years, we again observed a decline to- 
ward a minimum and a raise toward a new phase of strong 
emission, in a shorter period of about 3 years. 

As concerns, V/R ratio, peaks separation and radial 
velocity of central depression, we also observed a long 10 
years cycle out-of-phase with respect the EW and a more 
rapid variability occured during the last 3 years. As for the 
FWHM, we did not observe any appreciable change, at least 
at the resolving power of our data. 

Making use of some useful literature correlations be- 
tween FWHM, peaks separation and stellar projected rota- 
tional velocities, we estimated for the first time the vsini of 
(3 Cep Aa equal to ~ 230 km s _1 and the outer radius of the 
surrounding disk, rd~ 3.38 R». 
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APPENDIX A: JOURNAL OF OBSERVATIONS 



Table Al. Journal of observations. For each spectrum, ordered 
by heliocentric julian date (HJD) wc report: the equivalent width 
expressed in A, the ratio between V and R peaks, the peaks sep- 
arations and the radial velocity of the central depression. The 
last two quantities are all expressed in km s — 1 . The capital letter 
next to the HJD is a flag indicating the observatory: R - Rittcr 
Observatory, I - INT, E - Elodie and O - Catania Astrophysical 
Observatory. 



HJD 


EW 


V/R 


AV 


RVcd 


(2400000.+) 


(A) 




(km s- 1 ) 


(km s _1 ) 



49254 


.695 


R 


— 11. 


.5 


± 


1 


2 


1 


.06 


± 


0. 


10 


256 ± 13 


12 


1 


± 


1.4 


49261 


590 


R 


— 10. 


3 


± 


1 


3 


1 


.07 


± 


0. 


10 


250 ± 10 


18. 


.6 


± 


1.4 


49263 


.625 


R 


— 11. 


3 


± 


1 


1 


1 


.01 


± 


0. 


08 


255 ± 10 


16. 


.8 


± 


1.1 


49266 


.582 


R 


—9. 


9 


± 


3 


1 


1 


.08 


± 


0. 


10 


248 ± 13 


36 


.4 


± 


2.0 


49268 


.583 


R 


— 10. 


3 


± 


4 


,5 


1 


.04 


± 


0. 


10 


253 ± 12 


14 


2 


± 


2.0 


49271 


.583 


R 


-10. 


9 


± 


3 


1 


1 


.00 


± 


0. 


10 


250 ± 13 


33 


.5 


± 


2.3 


49319 


.551 


R 


— 10. 


1 


± 


4 


1 


1 


.06 


± 


0. 


14 


239 ± 16 


37. 





± 


3.1 


49333 


.625 


R 


—9 


.7 


± 





,8 


1 


.04 


± 


0. 


06 


247 ± 6 


17 


2 


± 


0.5 


49334 


.658 


K 


—9 


.4 


± 


1 


.4 


1 


.03 


± 





08 


248 ± 7 


20 


.8 


± 


0.8 


49579 


.878 


R 


—5 


9 


± 


1 


,7 


1 


.09 


± 





17 


249 ± 12 


24 


.1 


± 


2.5 


49582 


.721 


K 


—6 


.0 


zb 


2 


4 


1 


.05 


± 





13 


242 ± 11 


23 





± 


2.2 


49628 


.705 


R 


—6. 


7 


± 


1 


,5 


1 


.09 


± 


0. 


11 


247 ± 8 


26. 


.8 


± 


1.1 


50030 


.278 


E 


—4 


.8 


± 


1 


,7 


1 


.13 


± 


0. 


14 


250 ± 8 


39 


.8 


± 


1.2 


50710 


.563 


E 


-7 


.0 


± 





.9 


1 


.14 


± 


0. 


Id 


261 


± 


i I 


51 


,7 


± 


2.1 


50768 


.467 


1 


—5. 


9 


± 





.4 


1 


.20 


± 


0. 


15 


264 


± 


11 


42 


3 


± 


1.3 


50770 


.489 


I 


-5. 


.9 


± 


2 


,8 


1 


.16 


± 


0. 


19 


268 


± 


15 


52. 


.4 


± 


3.3 


51144 


,261 


E 


— 0. 


.6 


± 





.6 


1 


.79 


± 


0. 


75 


z62 


± 


26 


58 


.8 


± 


4.0 


51155 


.338 


I 


— 


.5 


± 





,7 


2 


.27 


± 


1 


04 


268 


± 


32 


76 


.9 


± 


4.0 


51325 


.682 


1 


— 


.6 


± 


1 


2 


1 


.66 


± 


1. 


.38 


281 


± 


28 


72 


.6 


± 


4.3 


52092 


.823 


R 


— 11. 


3 


± 





.8 


1 


.01 


± 


0. 


09 


260 


± 


14 


2 


.3 


± 


1.7 


52104 


.862 


R 


— 10. 


6 


± 





.6 


1 


.02 


± 


0. 


11 


261 


± 


14 





.9 


± 


3.2 


52155 


.713 


K 


—9 


9 


± 





,5 


1 


.03 


± 


0. 


09 


255 


± 


10 


2 


.8 


± 


1.6 


52157 


.795 


R 


—9. 


6 


± 


0. 


4 


1 


.05 


± 





14 


257 


± 


16 





.6 


± 


3.6 


52158 


.715 


R 


—9 


9 


± 


1 


1 


1 


.01 


± 


0. 


10 


256 


± 


12 


11. 


.0 


± 


1.9 


52159 


624 


K 


—9 


.4 


± 


3 


3 


1 


.02 


± 





11 


258 


± 


12 


22 


.3 


± 


2.0 


52167 


.710 


R 


—9 


7 


± 





5 


1 


.05 


± 





11 


257 


± 


11 


—2 


.4 


± 


1.9 


52220 


580 


R 


—9 


6 


± 


2 


4 





.99 


± 


0. 


11 


259 


± 


11 


16 


.9 


± 


1.9 


52223 


.565 


R 


—9 


6 


± 


1 


6 





.97 


± 


0. 


12 


241 


± 


15 


3 


.1 


± 


3.0 


52861 


.537 


E 


—9 


6 


± 


1 


2 


1 


.03 


± 





12 


255 


± 


15 


27 


1 


± 


2.0 


52862 


.436 


E 


—8. 


.6 


± 





.1 





.97 


± 





11 


262 


± 


14 


-6 


.3 


± 


2.0 


52864 


.634 


E 


-10. 


2 


± 


1 


.2 





.96 


± 


0. 


10 


258 


± 


14 


15. 


.9 


± 


1.8 


52871 


.460 


E 


-9 


.3 


± 





,8 





.95 


± 


0. 


10 


251 


± 


17 


11 


.7 


± 


2.0 


53240 


.765 


R 


-2. 


.2 


± 


4. 


.2 


1 


.06 


± 


0. 


36 


255 


± 


19 


39. 


.4 


± 


5.4 


53250 


.731 


R 


-1. 


3 


± 


5 


,1 


1 


.03 


± 


0. 


26 


240 


± 


18 


47 





± 


5.4 


53261 


.625 


R 


-2 


.0 


± 


2 


.4 


1 


.07 


± 


0. 


33 


255 


± 


16 


30 


6 


± 


3.8 


53264 


.599 


R 


-0. 


.8 


± 


5 


.0 


1 


.19 


± 


0. 


49 


261 


± 


19 


53 


6 


± 


5.3 


53271 


.686 


R 


-2. 


6 


± 


3 


,7 


1 


.25 


± 





36 


251 


± 


15 


d6 


,4 


± 


3.2 


53276 


.663 


R 


-2 


,2 


± 





,8 


1 


.10 


± 


0. 


2d 


256 


± 


14 


37. 


3 


± 


2.4 


53278 


.598 


R 


-2. 


.8 


± 





9 


1 


.31 


± 


0. 


3d 


243 


± 


14 


44 


.9 


± 


1.8 


53279 


.627 


R 


-2 


,5 


± 


2 


9 


1 


.26 


± 


0. 


3d 


238 


± 


16 


45. 


.5 


± 


3.2 


53281 


.697 


R 


-2 


9 


± 


2 


3 


1 


.08 


± 





18 


230 


± 


10 


30 


.8 


± 


2.1 


53282 


.588 


R 


-2. 


,1 


± 


1 


,4 


1 


.17 


± 





.42 


250 


± 


18 


25 


,5 


± 


5.4 


53284 


.703 


R 


-2 


9 


± 


1 


,1 


1 


.03 


± 


0. 


44 


240 


± 


28 


31 





± 


12. 


53286 


.550 


R 


-2. 


.5 


± 





,6 


1 


.10 


± 


0. 


35 


244 


± 


17 


36. 


,1 


± 


3.5 


53289 


,542 


R 


-2 


,5 


± 





,5 


1 


.08 


± 


0. 


31 


240 


± 


15 


29. 


.5 


± 


2.7 


53310 


.608 


R 


-1 


.5 


± 





9 


1 


.14 


± 


0. 


35 


244 


± 


15 


32 


.9 


± 


2.4 


53315 


.542 


R 


-2. 


3 


± 





3 


1 


.12 


± 





32 


235 


± 


15 


34 


,5 


± 


2.4 


53316 


.540 


R 


-1 


.5 


± 


1 


.0 


1 


.17 


± 


0. 


40 


246 


± 


17 


36 


.5 


± 


2.6 


53318 


.598 


R 


-2. 


.2 


± 


2 


.2 


1 


.09 


± 


0. 


34 


247 


± 


15 


24. 


.1 


± 


2.7 


53322 


.619 


R 


-2. 


3 


± 


2 


3 


1 


.25 


± 


0. 


48 


243 


± 


19 


25. 


4 


± 


3.3 


53323 


.620 


R 


-1 


,5 


± 


4. 


,5 





.92 


± 


0. 


37 


245 


± 


21 


-12 


.2 


± 


6.1 


53324 


.524 


R 


-2 


.5 


± 


1 


.5 


1 


.03 


± 


0. 


38 


252 


± 


20 


13. 





± 


5.0 


53336 


.500 


R 


-1. 


.4 


± 





.6 


1 


.17 


± 


0. 


40 


244 


± 


16 


37 


.4 


± 


2.6 
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Table Al - continued 



HJD 

(2400000.+) 




EW 

(A) 






V/R 




(km 


AV 

s-i) 


RVcd 

(km s- 1 ) 


53579.505 


O 


-0. 


.3 


± 


1 


.4 





.87 


± 


0. 


76 


227 


± 


35 


-29.3 


± 


11. 


53580.546 


O 


-0 


.2 


± 





.7 





.87 


± 


0. 


81 


245 


± 


30 


-38.1 


± 


7.0 


53581.559 


o 


2 


.2 


± 


1 


9 





.78 


± 


0. 


53 


248 


± 


29 


-35.4 


± 


4.2 


53635.525 


o 


2 


.2 


± 


1 


.2 





.98 


± 


0. 


.84 


261 


± 


34 








53672.415 





1 


,7 


± 


6 


.3 





.71 


± 


0. 


61 


264 


± 


48 








53673.457 


o 


-1. 


.5 


± 





.8 





.59 


± 





31 


266 


± 


36 








53674.475 








.8 


± 


4 


.7 





.84 


± 


0. 


54 


236 


± 


51 








53871.463 





-3 


.5 


± 





,7 





.96 


± 


0. 


26 


287 


± 


20 


-18.3 


± 


3.9 


53872.511 





-4. 


.6 


± 


2 


6 





.95 


± 


0. 


35 


298 


± 


24 


-14.1 


± 


4.6 


53873.522 


o 


-3 


.8 


± 


2 


.0 


1. 


.03 


± 


0. 


41 


288 


± 


26 


-5.9 


± 


5.5 


53933.850 


o 


-3 


6 


± 


5 








.92 


± 


0. 


27 


273 


± 


22 


-40.3 


± 


3.3 


54013.295 





-4. 


.6 


± 


5 


.0 





.97 


± 





23 


280 


± 


21 


-19.8 


± 


5.8 


54016.544 


o 


-5. 


3 


± 


5. 


.9 





.93 


± 


0. 


.20 


296 


± 


24 


-19.1 


± 


5.2 


54017.468 


o 


-4. 





± 


6 


.1 





.94 


± 


0. 


.21 


285 


± 


22 


-18.4 


± 


4.9 


54018.399 





-5. 


.5 


± 


4. 


.8 





.86 


± 





.18 


289 


± 


21 


-20.5 


± 


4.7 


54019.552 





-8 


.0 


± 


4 


.5 





.89 


± 





19 


297 


± 


18 


-23.1 


± 


3.5 


54021.544 





-5. 


.6 


± 


2 


.2 





.93 


± 


0. 


17 


275 


± 


17 


-30.4 


± 


2.9 


54022.398 


o 


-5. 


.2 


± 


5. 


.5 


1. 


.08 


± 


0. 


.22 


269 


± 


17 


16.9 


± 


3.5 


54041.432 





-4. 


6 


± 


3 


.0 





.96 


± 





.17 


274 


± 


16 


6.8 


± 


2.9 



